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identify and model molecular hydrogen absorption features longward of the Lyman 



> ■ edge at 912 A are described. Such features are commonly encountered in spectra ob- 

^ ■ tained by the Far Ultraviolet Spectroscopic Explorer and also in spectra obtained by the 

' Space Telescope Imaging Spectrograph, albeit less commonly. Individual templates are 

. calculated containing all the Lyman and Werner transitions originating from a single 

. rotational state {J") of the 0th vibrational level {v") of the ground electronic state. 

' Templates are provided with 0.01 A sampling for doppler parameters ranging from 

2 < 6 < 20 km and rotational states < J" < 15. Optical depth templates for 
excited vibrational states are also available for select doppler parameters. Each tem- 



^ ' plate is calculated for a fiducial column density of log[N(cm ^)] =21 and may be scaled 

^ ■ to any column less than this value without loss of accuracy. These templates will fa- 

cilitate the determination of the distribution of molecular hydrogen column density as 



a function of rotational level. The use of these templates will free the user from the 
computationally intensive task of calculating profiles for a large number of lines and 
\ allow concentration on line profile or curve-of-growth fitting to determine column den- 

sities and doppler parameters. The templates may be downloaded freely from the URL 
http:/ / www.pha.jhu.edu / ~stephan /h2ools2.html. 

Subject headings: ISM: abundances - ISM: molecules - line: identification - line: profiles 
- methods: analytical - molecular data - ultraviolet: ISM 



1. Introduction 

The launch of the Far Ultraviolet Spectroscopic Explorer FUSE has opened a new spectral 
window on the universe with a 905 - 1187 A bandpass (Sahnow et al. 2000). Within this window 
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spectra of continuum emitting objects almost universally exhibit absorption bands of molecular 
hydrogen, indicating this gas resides in the interstellar medium somewhere between the object and 
the observer. These bands arise from the photo-excitation of molecular hydrogen by the background 
continuum object, causing electrons in the ground electronic state X to be excited into either 
the higher B (Lyman band), or C^II^ (Werner band) states. The strength and density of the 
absorption lines are proportional to the distribution of molecular hydrogen column densities among 
the ground state rotational and vibrational (J", v") energy levels^ peculiar to the line-of-sight. 

For some observers the molecular hydrogen electron population distribution within the ro- 
vibrational states is a subject of direct interest. For others, molecular hydrogen absorption is a 
nuisance, potentially contaminating or even obscuring spectral features of interest. Fortunately 
there is a high degree of correlation between the absorption lines within different ro-vibrational 
bands {v'- 0) imposed by the uniformity of the energy level separations and the slow monotonic 
variation of oscillator strength as a function of upper vibrational state, a consequence of the Franck- 
Condon principle (Herzberg 1950). This correlation greatly aids the identification of molecular 
hydrogen features and provides a means for their removal from blended features, provided they can 
be reliably modeled in an unobscured portion of the spectrum. 

Here we will discuss the calculation and use of a simple set of optical depth templates from 
which observers can easily identify and model the molecular hydrogen features commonly observed 
in the FUSE spectral range. Each template contains all the Lyman and Werner absorption lines 
arising from the photo-excitation of electrons out of common ground state rotational J" and vi- 
brational v" = levels (usually) and into the v' states of the upper bands {v' <— 0), subject to 
the appropriate rotational state selection rules. A single column density N{v", J") controls the 
line strength and the templates are scaleable with this quantity. The optical depth variation with 
wavelength for each line is calculated using a Voigt profile to allow for the proper variation of 
line shape with increasing column density, reproducing the shapes endemic to the linear, flat, and 
square root portions of the curve-of-growth as controlled through the specification of the doppler 
parameter (b). The templates are given as optical depths functions, as opposed to transmission 
functions, to facilitate the inclusion of various atomic and molecular species or multiple velocity 
structures, simply by adding additional optical depths to the template function at the appropriate 
wavelength position. 

These templates can be used in conjunction with a variety of methods to determine the ap- 
propriate column densities and doppler parameters. One method is direct profile fitting, where the 
templates for various N{v", J") and b are used with a predetermined continuum to directly compare 
to the observed spectrum in a calculation. Another method uses the curve-of-growth technique, 
where equivalent width measurements are used to determine the appropriate column and doppler 



^When referring to the rotational and vibrational quantum numbers in a transition going from an upper electronic 
state (e.g. either B or C ^Ilu) to the ground state {X it is traditional to designate the upper state numbers 
with a single prime (e.g. J', v') and the lower state with a double prime (e.g. J", v"). 
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parameter from which the absorption Une transmission function can be quickly calculated. Once 
the columns, doppler parameters and continua are known in an unblended portion of the spectrum 
the degree to which they effect the other features, such as O VI or D I, maybe assessed. 

The template calculations incorporate the Abgrall et al. (1993a, b) databases of transition 
lifetimes, wave-numbers, rotational and vibrational quantum numbers for Lyman and Werner band 
transitions. The templates are available for free access on-line through http : //www.pha. jhu. edu/~stephan/h2ool 
Here we provide potential users with a ready reference to the calculation procedures, some examples 
of use, and discuss appropriate positional accuracy and column density limitations. We include for- 
mulae for the calculation of energy levels, population distributions, and line profiles in Appendices 
A and B. 

2. Energy Levels, Population Distribution, and Selection Rules 

The ground state rotational energy levels in the v" = vibrational state have temperature 
equivalent energies B^hc/k ^ 0, 170, 510, 1014, 1680, 2500 K for J = 5, while the 

vibrational energy levels are much more widely spaced with separations A G"hc/k = {G{v" + 1) — 
G{v")){hc/k) < 5983 K (see Appendix A). Boltzmann's law shows that the first few rotational 
states are low enough to be populated at the relatively cold temperatures associated with molecular 
gas in the interstellar medium, while excitation of the upper vibrational levels requires a more 
energetic environment. Absorption from several rotational levels is the norm, although the early 
Copernicus observations showed the typical distribution of column density in the higher rotational 
states is larger than expected from a single temperature Boltzmann distribution (Aannestad & 
Field 1973). The reasons for the deviations from a Boltzmann distribution are myriad and we 
refer to the reader to the comprehensive review by Shull & Beckwith (1982) who have discussed 
the various diagnostics for determining formation rates, collision rates and local radiation field 
intensity by using the rotational column density distribution of molecular hydrogen and the atomic 
hydrogen column. Absorption from excited vibrational states are rare. When present they indicate 
that the gas is affected by either a high temperature process such as shocks (McCandliss 2001) or 
a nonthermal processes such as fluorescent pumping (Meyer et al. 2001). 

Lyman bands have two branches R(J"), with J' = J" + I and P(J"), with J' = J" - 1. The 
selection rule that requires A J = it 1 results from the zero change in electronic angular momentum 
between the X^S^ and B^T,'^ states (Herzberg 1950). The Werner bands have an additional 
branch Q(J") with AJ = 0, allowed because the change in electronic angular momentum between 
the X'^T.+ and C^H^ states is one. In addition, the minimum J' for the state is 0, while 

the minimum J' for the C ^liu state is 1. This leads to an absence of P(0) in the Lyman bands and 
P(0), P(l) and Q(0) in the Werner bands. 

There are 377 absorption line transitions from Lyman (0 < v' < 20) and Werner (0 < v' < 5) 
bands arising from rotational levels, < J" < 6 in the ground vibrational level {v"= 0) of the 
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X^Tj'^ electronic state longward of the Lyman edge. Many more are possible if v" > becomes 
populated. Consequently the identification of the molecular hydrogen lines can be daunting. It is 
typically the first task in the post pipeline analysis of FUSE spectra. The identification process 
is eased considerably by locating the strongest Lyman lines usually R(0) for J" = and P(l), 
R(l) for J" = 1, and then proceeding to higher J" until the lines merge with the noise. Our 
computation of templates containing only those lines originating from a single rotational state 
allows for a direct implementation of this method. It greatly facilitates the determination of the 
N{v" , J") distribution, and allows for assessing the degree of contamination in coincident alignments 
of molecular hydrogen absorptions with lines from other diagnostically useful atoms or molecules. 

3. Optical Depth Template Calculation 

A total optical depth function is constructed for the J" rotational state of the v"= vibrational 
level from the sum of all the individual optical depth profiles ranging over all v' states in the Lyman 
{0 < v' < 18) and Werner (0 < v' < 5) line database given as, 

Tj"(A) = s:;r=»jr,(A). (1) 

The function rj//(A) is established on a common wavelength grid spanning the region of interest 
(790 - 1490 A in our case) with 0.01 A spacing. This total optical depth function is used to 
accumulate individual optical depth profiles. Each line profile is calculated using the Voigt function 
formulae in Appendix B, covering a range of wavelengths spanning 50 A to either side of the line 
center, so as to accommodate the broad Lorentzian wings that appear in the line profiles at high 
column density. The span of the common wavelength grid is arbitrary. It allows for the convenient 
inclusion of contributions from lines of absorbing species both inside and outside of the FUSE 
bandpass. An insertion point is determined to be the wavelength grid point closest to the line 
center and the individual optical depth profile is added to the total optical depth about that point 
and the process is repeated for all Lyman and Werner transitions arising from the same lower level 
rotational state. 

For the J" = states the template contains only the R(0) lines of the Lyman and Werner 
bands of {v' ^ 0). For J" = 1 there are R(l) and P(l) lines in the Lyman bands and R(l) and 
Q(l) lines in the Werner bands. Above J" = 1 there are R(J") and P(J") lines in the Lyman 
bands and R(J"), Q(J") and P(J") lines in the Werner bands. A representative set of templates 
for < J" < 4 with a doppler parameter of 5 km s~^ with log columns of 21 and 20 for J" = 0,1, 
19 to 18 for J" = 2, 17 for J" = 3, and 16 for J" = 4. are displayed in Figure 1. 

4. Line Position, Equivalent Width and Column Density Limits 

The location of each line within the larger wavelength grid is accurate to within the grid 
step or 3 km s^^. This method of insertion is computationally expedient compared to the more 
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accurate but much slower method of interpolating the centered profile onto the 0.01 A grid. Our 
grid spacing is on the order of the largest differences between observed and calculated wave-numbers 
in the Abgrall et al. (1993a, b) databases (~ 1 cm~^). The resolution of the FUSE spectrograph is 
R ~ 20000 or 15 km s~^ (Sahnow et al. 2000), so the grid oversamples by a factor of 5. 

The minimum grid step of 0.01 A sets a lower limit on the useful doppler parameter range. The 
core of a line with a doppler parameter of 1 km s~^ is not sampled with sufficient resolution to yield 
an accurate measure of optical depth with increasing distance from the line core. Consequently, the 
numerically integrated equivalent width from the 1 km s~^ template is higher than the theoretical 
value in the linear portion of the curve-of-growth and lower in the flat portion of the curve-of-growth. 
Although available on the website, the 1 km s~^ templates should not be used for quantitative work. 
This problem is quantified in Figure 2 where we show curves-of-growth (in red) for a range of NfX, 
and a descrete set of doppler parameters (b=l, 2, 4, 8, and 16 km s~^). Overplotted are the 
equivalent widths derived from direct numerical integration of the Lyman lines in the J=0 and J=l 
templates over the same range of doppler parameters. The problem with the equivalent widths 
calculated from the lines contained in the b = 1 km s~^ template is evident. There is a slight 
(< 7%) residual from the theoretical curves in b = 2 km s^^ derived equivalent widths. 

An examination of the current literature shows that very few lines-of-sight have been found 
with molecular hydrogen determined doppler parameters < 2 km s~^. For example, in a study of 
the Large and Small Magellanic clouds Tumlinson et al. (2002) find only two lines-of-sight (out 
of 37 determinations) that have molecular hydrogen doppler velocities of < 2 km s~^, although 
several have lower bounds that dip below b = 2 km s~^. Another study by Meyer et al. (2001) 
finds a b = 1.8 (±0.1) km s~^ for the highly excited molecular hydrogen in the line-of-sight toward 
HD 37903, the central star of the reflection nebula NGC 2023. These instances notwithstanding, 
the current set of templates should be adequate for most work. We note that lines-of-sight with 
b = 2 km s~^ will have instrumental line profiles, except where column densities are high enough 
to produce damping wings. The curve-of-growth technique will be prefered to line profile fitting for 
determining column densities and doppler parameters in these instances, as noted in the following 
section. The templates are not used in the curve-of-growth analysis except to represent the solution 
after it has been found. The regime where absorption lines will have instrumental profiles in the 
FUSE spectrographs is roughly identified in Figure 2 as the region below the solid horizontal line 
drawn at logW/X = 1/R. 

Nevertheless there are many cases where small doppler velocities have been measured in species 
that correlate well with molecular hydrogen, as in the case of CH (e.g. Andersson, Wannier, & 
Crawford (2002) and Federman (1982)). In order to support such investigations in the future, 
templates will be developed with finer wavelength and doppler velocity grids for exploring the low 
b regime. Higher resolution grids will also find use in studies requiring closely spaced multiple 
velocity components with small doppler velocities. However, the current set of grids are large (~ 
8Mbyte uncompressed for a single b up to J" = 15) and the high resolution grids will be at least 
twice as large. 
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The column densities have no lower limits, they may be scaled to as low a column as required 
without loss of accuracy. An upper limit results from the finite wavelength span of ±50 A for the 
individual profile calculations. Above logN{v", J") = 21 (the fiducial column density for which 
each optical depth template is calculated) the individual line profiles suffer a progressively severe 
truncation of the Lorentz wings at ±50 A from line center, which manifests itself as a series of 
discontinuities in the summed optical depth grid. The log A^(t;", J") < 21 is a "fuzzy limit" as the 
seriousness of the truncation depends upon the oscillator strength. This effect is independent of 
doppler parameter because at logA^(f", J") > 21 all lines are fully damped for b < 20 km s~^. 

The effect is illustrated in Figure 3 with the lines of minimum and maximum oscillator strength 
in the R(0) transitions of the Lyman bands (0 - 0) and (7-0), where we show the half profiles 
of these two lines over a 50 A span. We see that for a column density of log A^(0, 0) = 21 the line 
profile has returned from black at the line core to 0.9995 and 0.9990 at 50 A from line center for 
(0 - 0) and (7 - 0), respectively. For a column density of logA^(0, 0) = 22 the lines only return to 
0.995 and 0.990, while for logiV(0, 0) = 23 the return is 0.95 and 0.90, respectively. It is difficult to 
estimate quantitatively the combined effects of profile truncations since the error it generates is a 
function of wavelength. The general effect is to cause a progressive underestimation of the combined 
blackness of the rotational template for column densities log N(y", J") > 21. Quantitative work 
at columns much greater than \ogN{v", J") > 21 using these templates should be avoided. This 
should not be too much of a problem as few sight lines with logN{v", J") > 21 have been found. 
For instance, in an exploration of 23 sight lines with Ay ^ 1 Rachford et al. (2002) found only 1 
case (HD 154368) where the formal value of logiV(J" = 0) was 21.04 ± 0.05. 



5. Discussion and Examples 

The templates are stored in files with names that encode the calculation parameters. For 
example, tauh2n21b2j0-15v0.dat, indicates a column density of log iV(cm~^) = 21, a doppler width 
of 2 km s~^, a rotational state range of of < J" < 15, and a ground vibrational state of v" = 0. 
The files contain 17 arrays of unformatted double precision numbers with 59000 elements. The 
first array is the wavelength variable. The next 16 arrays are the J" =0, 1, • • • 15 optical depth 
functions. Table 2 gives a schematic of how the data arrays are placed in the file. A simple IDL 
program is provided on the website to read in the data and form an example transmission function. 
There are extra empty elements longward of the FUSE long wavelength cutoff at 1189 A in part 
because there are a few high rotational lines that exist beyond this cutoff. However, the main reason 
is to allow for the inclusion of spectra for ground state vibrational levels v" > 0, for the purpose 
of investigating molecular hydrogen excitation processes in energetic astrophysical environments. 
Rotational templates for v" > states will be included on the H200IS website in the near future. 

Absorption spectra (the transmission functions) are obtained by taking e"^^. All the optical 
depth templates are tabulated as negative numbers. The spectra for column densities other than 
10^^ particles cm~^ can be obtained by dividing the optical depth by 10^^ and multiplying by 
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the required column density value. For example, the absorption spectrum for 10 is given by 
grA/iooo These spectra should be convolved with an instrumental line spread function to simulate 
instrumental broadening before comparing to an object spectrum. For use with FUSE spectra, 
convolution with a gaussian width with a FWHM on the order R > 15000 is typical. For high 
precision work it may be necessary to employ a wavelength dependent convolution kernel, as the 
resolution of the various spectral channels is wavelength dependent at the 10% level (Sahnow et al. 
2000). 

Once convolved, an overlay of the J" = transmission function ( multiplied by an appropriate 
continuum) on an object spectrum will allow an assessment of the need for zero-point offsets to 
the wavelength registration. The J" = spectrum contains the bluest lines in any given Lyman 
band {v' - 0), so it should line up with the blue side of the blended R(0) and R(l) lines. For the 
highly damped lines encountered at high column density it may be necessary to include higher J" 
states to assess zero-point offsets. Adjust the column density of the absorption spectrum until a 
good fit is obtained. Once satisfied add the optical depth template from the next highest rotational 
state, take the exponent, adjust the columns and wavelength zero-point as necessary, and reassess 
the fit. Repeat until the inclusion of a higher level J" no longer improves the fit. The process 
may be quantified by computing the statistic over wavelength intervals selected to coincide with 
unblended absorption features in the object spectrum. Repeat for different column densities and 
doppler parameters and then find the minimum x^. 

It it is often expedient to perform the initial fitting operation by eye. Interactive assessment 
can be used to determine the columns to within a few dex and the doppler parameters within a few 
km s~^, provided complicated velocity structures on the line of sight are absent and the column 
density is high enough. These approximate determinations can serve as the starting point for more 
precise work and minimize the time devoted to minimization by providing good guesses. In 
Figure 4 we show how the higher vibrational transitions can be used to define the column and the 
lower vibrational transitions the doppler parameter when the absorptions are making the transition 
from the flat portion of the curve-of-growth to the square root portion where the lines profiles 
become fully damped. At log N{J" = 0) = 19 we see little difference in the higher vibrational line 
profiles for the two doppler parameters, while the lower vibrational lines do show a difference. At 
low doppler parameter the lines are still saturated, having box like profiles, while damping wings 
have developed in the higher doppler parameter line. 

In the low column transition region (logN(J" = 0) ~ 15) between the linear and flat portion 
of the curve-of-growth, discerning the tradeoff in column and doppler parameter is less apparent 
to the eye. In this regime the lines are unresolved (by FUSE) and will have proflles close to the 
instrumental profile except for the highest doppler parameters. In these cases the curve-of-growth 
method for determination of column is generally more useful than profile fitting. In between these 
two regimes on the flat (slowly varying) portion of the curve-of-growth, it is relatively easy to 
determine the doppler parameter (b) from the box like profiles, but the column is ill-defined. Here 
profile fitting is the most accurate method for column determination. 
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In Figure 5 we show a curve-of-growth example. The molecular hydrogen column densities 
and a doppler parameter have been derived from a curve-of-growth analysis of all the lines in the 
rotational states < J" < 5 between 1040 ~ 1098 A. The bottom right panel for the (4-0) Lyman 
band fit shows an overplot of the transmission function computed by applying the derived column 
densities and doppler parameter. A curve-of-growth analysis was then applied to eight D I lines 
starting at Lyman-7. An absorption line model containing lines from H I, D I, O I, and molecular 
hydrogen has been overplotted in each panel. 

In Figure 6 we show an example of a fit to the (4-0) Lyman band in the lower panel, which 
was fine-tuned by taking into account the fit to the (0-0) Lyman band. The rotational templates 
were used with an appropriate continuum to calculate for selected portions of the line profiles. 
In the upper panel of we illustrate how well the II2 model from the (4-0) band fit reproduces the 
H2 structures in the (5-0) band. The template model can be used to constrain the location and 
degree of blending with lines of O VI, C II, O I and Ar I. Several unidentified lines are also evident, 
some of which are higher rotational level II2 lines. 

These are just two examples illustrating how the use of these templates can expedite the 
identification, computation and deblending of molecular hydrogen profiles from far-UV spectra. The 
use of a common 0.01 A wavelength grid facilitates the inclusion of profiles from other atomic and 
molecular species. In addition, any number of multiple velocity components can easily be coadded 
and scaled, with the appropriate doppler dilatation or contraction applied through the interpolation 
of the optical depth template onto a shifted wavelength grid. Such interpolated optical depth 
templates can be employed to search for molecular hydrogen in intervening intergalactic absorption 
systems that appear in the spectra of quasars and galaxies at high redshift. The templates can 
find use in any study that requires an accounting for molecular hydrogen blending. Examples 
include studies seeking to determine the velocity distribution and abundance of O VI in the galaxy 
and halo, of D I in the local and nearby ISM, or of any other interstellar atomic and molecular 
species that may blend with molecular hydrogen features. The templates also provide a means 
to determine the excitation state of molecular hydrogen so as to directly test the predictions of 
various excitation mechanisms, such as fluoresence (radiative pumping), formation pumping, or 
shock excitation. These templates can easily be used with \^ minimization codes to assess the 
goodness-of-fit as functions of column density and doppler parameter. They provide a means 
for the user to concentrate on the fitting procedure without having to take on the computationally 
intensive task of calculating line profiles for a large number of absorbers. Such goodness-of-fit codes 
are under development and testing by many groups (c.f. Owens. f discussed in Lemoine et al. 2002 
and Rachford et al. 2002), and some of these codes are available on-line^. The author will entertain 
inquiries into the use of beta versions that he is developing as time allows. 

I would like to acknowledge David Neufeld for the original suggestion to calculate templates for 



^http: //origins. colorado.edu/~tumlinso/h2/sw/sw. html 



- 9 - 



the individual rotational states and William Blair for organizing a series of tutorials on molecular 
hydrogen that provided the impetus for this work. I would also like to acknowledge Brian Wolven 
for providing an electronic version of the Abgrall et al. ( 1993a, b) databases and Ravi Sankrit for 
commenting on initial drafts of this paper. The calculations described here were performed in IDL 
on a Unix workstation. Some users, working on non-unix machines, have reported problems with 
byte swapping when downloading the templates and I would like to thank Pat Romano for kindly 
providing a set of byte swapped templates that were developed for use on linux machines. Finally, 
simulating discussion provided by various members of the FUSE science and operations team is 
greatly appreciated. This work has been supported by NASA grants NAG5-5315, NAG5-10403, 
and NAG5-11456 to the Johns Hopkins University. 



A. Term Values and Thermal Population Distribution 

Following Herzberg (1950), the ro- vibration energy levels (in wave-number units of cm~^) of 
an anharmonic oscillator is a sum of electronic, vibrational and rotational term values, 

T{v, J) = Te + Giv) + F{v, J) (Al) 

where Tg is a constant for the given electronic state. G{v) and F{v, J) are given as expansions in 

{v + I) and J(J + 1) respectively with V = 0,1 ■■■ Vmax, J = Jmin, Jmin + 1 • • • Jmax, 

Giv) = tOeiv + ^) - LJeXeiv + + LOeVeiv + + ■ ■ ■ , (A2) 

F{v,J) = B^J{J+l)-D^{J{J + l)f + ■■■, (A3) 
and By and are also given as expansions in (f + ^), 

By = Be- ae{v + ]^) + --- , (A4) 

Dy = De + f3e{v + ]^) + ---. (A5) 

We reproduce in Table 1 the most significant term value constants for the X ^ , B ^ , and 
C^Iiu states collected by Huber &: Herzberg (1979). They have also compiled a comprehensive 
list of the higher order constants of these expansions. 

The X^Tj^ energy levels may be used in Boltzmann's law to determine the distribution of 
electrons among the ro-vibrational levels of the ground state under the assumption of thermal 
equilibrium. Following Draine & Bertoldi (1996) we write, 

N{v, J)/g.j = N expf^^^^^^^] (A6) 

where is a renormalization factor determined by first calculating the matrix elements N(v,J) 
with N = 1 and then defining N = N [v , J) for the renormalization. gj is the degeneracy 

of the states where for even J (para-H2) gj = 2 J + 1 and for odd J (ortho-H2) gj = 3(2J + 1) 
under the assumption that the ortho to para ratio is 3. 
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B. Absorption Profiles 

The formulae reproduced here are quite general and may be used in the calculation of any 
absorption line profile, atomic or molecular. The shape and strength of an absorption line depends 
on 5 independent constants, the oscillator strength iv'v"j'j"-, the wavelength Xv'v"j'j", the total 
lifetime of the upper state r^,/^,//j/j//, the doppler velocity b = ^J{2kT/mH2Y + '^fur-fe' ^^e 
column density of the lower state N(v",J"). The first 3 variables are set by quantum mechanics 
while the latter 2 are the quantities of interest along the line-of-sight. At first glance it would appear 
that we need to calculate separate templates for a grid of column densities and doppler parameters. 
However, this is not necessary because the cross-section profile is independent of column density 
so the optical depth scales linearly with column density. Consequently a series of optical depth 
templates, calculated with the same doppler velocity, can be directly scaled with column density, 
reducing our need for an optical depth series where the doppler parameter alone is varied. 

Following Cartwright & Drapatz (1970), the normalized line profile as a function of wavelength 
(A) for the absorbing transition of an electron in level v'J' <— v"J" (which we will designate with 
an index i) maybe expressed as: 





= exp(-Tj(A)), 


(Bl) 




= N{v",J")a,{X), 
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Here H{a, u) is the Voigt function, a convolution of the a Maxwellian velocity distribution with 
a Lorentz profile (Rybicki & Lightman 1979). The dimensionless frequency is n = and 
a = ^ ^ is the damping parameter. The continuous, discrete and doppler frequencies are 
u = c/A, Vi = c/Xi, and Vd = {b/cjvi (c is the speed of light, e is the electron charge in esu, 
rUe is the electron mass in grams, and b is defined above). The total lifetime of the upper state is 
a sum over all possible vibrational transitions out of the upper v' state and into the lower states 
< v" < 14 such that Tj = Y^jf, nA' where A- = ^Jj,t} ^ fi- The line core cross-section is 

OiiXi) as the Voigt function has a value of i?(0, 0) = 1 at line center (when a « 1 as is usually 
the case). (7j(A) and Ti(A) are the cross-section and optical depth profiles, respectively. 

In Figure 2 the curves-of-growth shown in red result from the numeric integration of 1 — 
exp (— rj(A)) over the complete wavelength interval of the profile, for a range of column densities, 
and descrete set of doppler parameters. The (well known) formulae for the linear, flat, and square 
root portion of the curves are overplotted (in black) to directly show the regions of validity for the 
approximations . 



(B6) 



c mech 



me C'^ 



These formulae are discussed in detail in (Spitzer 1978). It can be seen in Figure 2 that the 
formula for the flat portion of the curve-of-growth is low by approximately 6 - 8% with respect to the 
numeric calculation. An examination of (Miinch 1968) shows this descrepancy can be accounted for 
by keeping the second order term in the asymptotic expansion of ^ for large tq = ( ^^'^ ) ^ 1- 
The correction factor to ^ is then (1 + .2886/ In tq) and setting tq ~ 50 yields a 7% correction. 
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Fig. 1. — The transmission functions as derived for all the Lyman and Werner absorption lines 
originating from J" = 0,1,2,3,4 longward of 900 A, top to bottom respectively. The vibrational 
transitions [v' - 0) are labeled above the lines with the Werner band label above the Lyman. The 
log columns are 21 and 20 for J" = 0,1, 19 to 18 for J" = 2, 17 for J" = 3, and 16 for J" = 4. These 
columns were chosen to be representative of the monotonically decreasing distribrutions commonly 
derived from far UV spectra. 

Fig. 2. — Curve-of-growth diagram illustrating the regions of validity for each of the equivalent 
width approximation formulae. The numeric solution is in red, and the approximations are in black. 
Overplotted are the equivalent widths derived from direct numerical integration of the Lyman lines 
in the J=0 and J=l templates for b = 1, 2, 4, 8, and 16 km s"^ showing the problem with the 
b = 1 km s~^ template. The solid black horizontal line at the middle of the graph indicates roughly 
the boundary between resolved and unresolved (instrumentally broaded) line profiles, assuming 
the FUSE resolution is i? = 20000. The dashed lines at the bottom of the graph indicate the 3cr 
equivalent width detection limits for continuum signal-to-noise ratios of 10 and 30 per resolution 
element. The placement of these lines was determined with the formula W = 3{N/S){\/R). 

Fig. 3. — On the right is a blowup of the Lyman R(0)(7-0) line as inserted into the common 
wavelength grid showing the truncation edge at -50 A from line center. The log scalings of 21, 22 
and 23 are overplotted. On the left is the same for Lyman R(0)(0-0), which has a smaller oscillator 
strength. The truncation edge is not as deep in the line on the right. 

Fig. 4. — The variation of line profile for a fixed column of log A = 19 in the R(0) (3-0), (2-0), 
(1-0) and (0-0) Lyman lines with two different doppler parameter (b) of 5 km s~^ and 10 km s~^ 
are overplotted. The shortward lines are very nearly the same, having become nearly fully damped 
they differ only in the core, while the (0-0) line for b = 5 km is still on the flat portion of the 
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curve-of-growth. This illustrates how the higher vibrational transitions can be used to define the 
column and the lower vibrational transitions the doppler parameter b when the absorptions are 
making the transition from the flat portion of the curve-of-growth to the square root portion where 
the lines profiles become fully damped. 

Fig. 5. — Curve-of-growth fitting example. The H2 column densities and a doppler parameter have 
been derived from a curve-of-growth analysis. The fit for the (4-0) Lyman band fit is displayed in 
bottom right panel. An absorption line model containing lines from H I, D I, O I, and H2 has been 
overplotted in each panel. The D I and O I models were also derived from curve-of-growth analysis 
while the H I model was "x-by-eye" from the Lyman-7 profile, which exhibits a damping wing. 

Fig. 6. — lins profile fitting example. The rotational templates were used with an appropriate 
continuum to calculate for selected portions of the line profiles of the (4-0) Lyman band. The 
upper panel shows how well the H2 model from the (4-0) band fit reproduces the H2 structures in 
the (5-0) band and constrains the location and degree of blending with lines of O VI, C II, O I and 
Ar I. 



Table 1. Molecular Hydrogen Term Value Constants 



State 






UJeXe 


Be 


Oe 


De 




100089.8 
91700.0 




2443.77 
1358.09 
4401.21 


69.524 
20.888 
121.33 


31.362 
20.015 
60.853 


1.664 
1.1845 
3.062 


2.23x10-2 
1.625x10-2 
4.71x10-2 
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Table 2. Template File Arrays 
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.P(2)(18- 0) 
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Flux (xlO"" ergs cm'^ s"' A"') 



Flux (xlO"" ergs cm"^ s"' A"') 
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o 
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Q(5)(2- 0) 
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